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Abstract 



We study how the solution of the two-dimensional Dirichlet boundary problem 
for smooth simply connected domains depends upon variations of the data of the 
problem. We show that the Hadamard formula for the variation of the Dirichlet 
Green function under deformations of the domain reveals an integrable structure. 
The independent variables corresponding to the infinite set of commuting flows are 
identified with harmonic moments of the domain. The solution to the Dirichlet 
boundary problem is expressed through the tau-function of the dispersionless Toda 
hierarchy. We also discuss a degenerate case of the Dirichlet problem on the plane 
with a gap. In this case the tau-function is identical to the partition function of 
the planar large N limit of the Hermitean one-matrix model. 
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1 Introduction 



The subject of the Dirichlet boundary problem in two dimensions |1| is a harmonic 
function in a domain of the complex plane bounded by a closed curve with a given value 
on the boundary and continuous up to the boundary. The question we address in this 
paper is how the harmonic function in the bulk varies under a small deformation of the 
the shape of the domain. 

Remarkably, this standard problem of complex analysis possesses an integrable struc- 
ture [§, m which we intend to clarify further in this paper. It is described by a particular 
solution of an integrable hierarchy of partial differential equations known in the litera- 
ture as dispersionless Toda (dToda) hierarchy. Moreover, related integrable hierarchies 
arise in the context of 2D topological theories and just the same solution to the dToda 
hierarchy emerges in the study of 2D quantum gravity 0, || (we do not elaborate these 
relations in this paper). 

Let D be a simply connected domain in the complex plane bounded by a smooth 
simple curve 7. The Dirichlet problem is to find a harmonic function u(z) in D such that 
it is continuous up to the boundary and equals a given function u (z) on the boundary. 
The problem has a unique solution written in terms of the Green function G(z\,Z2) of 
the Dirichlet boundary problem: 



u{z) = -^<fu (Od n G{z,0\dt\, (1.1) 

Z7T J-y 



where d n is the normal derivative on the boundary with respect to the second variable, 
and the normal vector n always looks inside the domain, where the Dirichlet problem is 

posed. Equivalently, the solution is represented as u{z) — — <£ uo(^)d^G(z,^)d^, where 

dD is understood as 7 run anticlockwise with respect to the domain. 

The main object to study is, therefore, the Dirichlet Green function. It is uniquely 
determined by the following properties Jl|: 

(Gl) The function G(zi,Z2) — log^i — z 2 \ is symmetric, bounded and harmonic every- 
where in D in both arguments; 

(G2) G(zi, Z2) = if any one of the variables belongs to the boundary. 

The definition implies that G(zi, z 2 ) is real and negative in D. The Green function can 
be written explicitly through a conformal map of the domain D onto some "reference" 
domain for which the Green function is known. A convenient choice is the unit disk. Let 
f(z) be any bijective conformal map of D onto the unit disk (or its complement), then 



G(z l ,z 2 ) = log 



f(zi) ~ f(z 2 ) 



f( Zl )f(z 2 ) - 1 



X2) 



where bar means complex conjugation. Such a map exists by virtue of the Riemann 
mapping theorem 0. 

It thus suffices to study variations of the conformal map f(z) under deformations of 
the boundary. This problem was discussed in @, [§, where it was shown that evolution 
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of the conformal map under changing harmonic moments of the domain is given by the 
dToda integrable hierarchy^. The study of the Dirichlet problem approaches this subject 
from another angle. 

Our starting point is the Hadamard variational formula J/J. It gives variation of the 
Green function under small deformations of the domain in terms of the Green function 
itself: 

8G( Zl ,z 2 ) = ±- <f d n G(z 1 ^)d n G(z 2 ,0Sh(0m (1.3) 

Here Sh(^) is the thickness between the curve 7 and the deformed curve, counted along 
the normal vector at the point (67. We show that already this remarkable formula 
reflects all integrable properties of the Dirichlet problem. 

A smooth closed curve 7 (for simplicity, we may assume it to be analytic in order 
to have an easy sufficient justification of some arguments below) divides the complex 
plane into two parts having the common boundary: a compact interior domain D int , 
and an exterior domain D ext containing 00. Correspondingly, one recognizes interior and 
exterior Dirichlet problems. The main contents of the paper is common for both of them. 
To stress this, we try to keep the notation uniform calling the domain simply D. We 
will show that (logarithms of) the tau-functions, introduced in [0, [| and further studied 
in ||, for the interior and exterior problems are related to each other by a Legendre 
transform. 

The exterior Dirichlet problem makes sense when the interior domain degenerates into 
a segment (a plane with a gap). We will show that in this case a deformation problem is 
described by the dispersionless limit of the Toda chain hierarchy and discuss its relation 
to the planar limit of the Hermitean matrix model. 



2 Deformations of the boundary 

Let Dbea simply-connected domain in the extended complex plane bounded by a smooth 
simple curve 7. Consider a basis ipk{ z ), k > 1, of holomorphic functions in D such that 
ipkizo) — for some point z G D. We call z the normalization point. The basis is 
assumed to be fixed and independent of the domain. For example, in case of the interior 
problem one may assume, without loss of generality, that the origin is in D and set 
zq = 0, ipk{z) = z k /k, while a natural choice for the exterior problem is zq = 00 and 
ipk{z) = z~ k /k. Throughout the paper, these bases for interior and exterior problems are 
refered to as natural ones. 

Let tk be moments of the domain D defined with respect to the basis ipk'- 

t k = k- [ Mz)d 2 z, k = 1,2,... (2.1) 

7T JD 

where k — ± for the interior (exterior) problem. We also assume that the functions ipk 
for domains containing 00 are integrable, or the integrals are properly regularized (see 
below). Besides, we denote by t the area (divided by tt) of the domain D in case of the 

1 A relation between conformal maps (of slit domains) and special solutions to some integrable equa- 
tions of hydrodynamic type was earlier observed by Gibbons and Tsarev B . 
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interior problem and that of the complementary (compact) domain in case of the exterior 
problem: 

(lr 

— d 2 z for compact domains 

7T JD 

to = I 

— / d 2 z for non-compact domains 

7T JC\D 

Let us note that the moments (except for to) are in general complex. We call the quantities 
tk, tk and to harmonic moments of the domain D. The Stokes formula represents the 
harmonic moments as contour integrals 



f k = 77— f ^k(z)zdz , k = 0,1,2,... 

ITXI J~i 



(where it is set ipo(z) = 1) providing, in particular, a regularization of possibly divergent 
integrals ( |2.1| ) in case of the exterior problem. Throughout the paper the contour in <f 7 
is run in anticlockwise direction both for interior and exterior problems. 

The basic fact of the theory of deformations of closed analytic curves is that the (in 
general complex) moments t^ supplemented by the real variable t form a set of local 
coordinates in the space of smooth closed curves || (see also fill). This means that 
under any small deformation of the domain the set {t ,ti,...} is subject to a small 
change and vice versa. More precisely, let 7(2) be a family of curves such that d t t^ = 
in some neighborhood of t = 0, then all the curves 7(i) coincide with 7 = 7(0) in this 
neighborhood. 

The family of differential operators 

V(s) = d t0 + J2(Mz)dt k +M*)dt k ) (2-2) 

k>l 

span the complexified tangent space to the space of curves. They are invariant under 
change of variables in the following sense: let tk be harmonic moments defined with 
respect to another basis, of holomorphic functions in D; then V(z) = V(z). Note 
that V(zo) = dt since ipk{zo) — 0. The operator V(z) has a clear geometrical meaning 
described below. 

Let us consider a special deformation of the domain obtained by adding to it an 
infinitesimal smooth bump (of an arbitrary form) with area e located at the point (67. 
Our convention is that e > if the bump looks outside the domain in which the Dirichlet 
problem is posed, as is shown in Fig. 0. 

Let A be any functional of a domain that depends on the harmonic moments only. 
The variation of such a functional in the leading order in e, is given by 

%)A = /^V(£M, £e 7 , (2-3) 

Indeed, combining 5 A = d to A5t + '^2(d tk A5tk + di k A5tk) and 

k>l 

If e 

St h = K / 1p k (z)d 2 Z = K-1pk(0 
VT J bump 7T 
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Figure 1: Action of the operator V(£) in the case of interior D^t (left) and exterior D ext 
(right) domains. In our convention bump always looks outside the (interior or exterior) 
domain. 



we obtain ( |2.3| ). So, the result of the action of the operator V(£) with £ 6 7 on A is 
proportional to the variation of the functional under attaching a bump at the point £. 
To put it differently, we can say that the boundary value of the function V(z)A is given 
by the l.h.s. of (|2.3| ). For functionals A such that the series V(z)A converges everywhere 
in D up to the boundary, this remark gives a usable method to find the function V(z)A 
everywhere in the domain. This function is harmonic in D with the boundary value 
determined from ( 2T3 ). It is given by ( |1 . 1| ) : 



V(z)A = «-L j \dt\d n G(z, 0^(0 A (2.4) 



This gives the result of the action of the operator V(z), when the argument is anywhere 
in D. 

For example, given any regular function / in a domain containing the interior domain 
D, set A f = J D f(z)d 2 z. We have: 5 e{0 A f = -V(£)A f = e/(f), £ e 7. If the function / 
is harmonic in D, then V(z)Af = 7if(z) for any z G D. 

The subject of the deformation theory of the Dirichlet boundary problem is to com- 
pute V(z)G(zi, Z2) through the conformal map or the Green function of the original 
domain. In the next section, we do this using the Hadamard variational formula. 



3 Hadamard variational formula and dispersionless 
integrable hierarchy 

3.1 The Hadamard integrability condition 

Variation of the Green function under small deformations of the domain is known due to 
Hadamard |7|], see eq. ( |1.3[ ). Being specified to the particular case of attaching a bump 
of the area e, it reads: 

8 m G( Zl ,z 2 ) = - ^d n G( Zl ,0d n G(z 2 ,0, f e 7- (3.1) 
To find how the Green function changes under a variation of the harmonic moments, 



we use (2.4) to employ the harmonic continuation procedure explained in the previous 
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section. The harmonic function in D with a boundary value given by the Hadamard 
formula is 

V(z 3 )G( Zl ,z 2 ) = k^- <f d n G(z 1 ,Od n G(z 2 ,Od n G(z 3 ,Om (3.2) 



It is obvious from the r.h.s. of (|3.2| ) that the result of the action of the operator V(z) 
on the Green function is harmonic and symmetric in all three arguments, i.e., 

V(z 3 )G( Zl , z 2 ) = V( Zl )G(z 2 , z 3 ). (3.3) 

This is our basic relation. It has the form of integrability condition. In the rest of the 
paper we will draw consequences of this symmetry and underlying algebraic structures. 
Note also that despite the Green function vanishes on the boundary, its derivative (the 
l.h.s. of eq. Q3.3p) with respect to the deformation of the domain does not. 

The basic equation (|3.3| ) is a compressed form of an integrable hierarchy. To unfold 
it, let us separate holomorphic and antiholomorphic parts of this equation. 

Let E be the exterior to the unit disk. Given a point a G D, consider a bijective 
conformal map f a :D^E such that f a {a) = oo. The Dirichlet Green function then is 

G(a,z) = -\og\f a (z)\. (3.4) 

Under a proper normalization of the map the integrability condition ( |3.3| ) becomes holo- 
morphic: 

V(&)]og/«Cz) = V(a)log/ 6 (*) (3.5) 

for all a,b,z G D. 

The following normalization will be convenient: the overall phase is chozen to be 
arg/ a (zq) = tt — arg (z — a) if a ^ z , where z G D is the normalization point. If a — z 
we set lim y(z — z ) 2 f'(z)^j to be real and negative. Under these conditions 

Ja{z) = -== -— — — (3.6) 

\(a-z )f(a)J f{a)-f(z) 

(for a, Zq 7^ oo). In the vicinity of the point a (a ^ oo) 

f a {z) = c **mo) Jf_ (i + Y /Pk (a)(z-a) k ), (3.7) 



z — a 



k>l 



where the real constant r a is called conformal radius of the domain fl]J with respect to 



the point a, and to is a phase determined from the normalization condition. In particular, 
they read u(z , z ) = 0. Similarly, the map f a (z) can be defined in the case when either 
a or Z q lies at infinity. 

To verify ( |3.5| ), we note that the holomorphic function V(6) log/ a (^) — V(a) log fb(z) 
is also antiholomorphic (in z) by virtue of (|3.3| ), and thus must be a constant. Setting 
z = zq we find that the latter is zero: 

V(6) log|/ a (z )|-V(a)log|A(z )| + 

+ zV(fo) arg/ a (^ ) - iV(a) &ig f b (z ) = 

(the first line vanishes due to ( |3.3|) , the second one vanishes because the normalization 
does not depend on the shape of the domain) . 
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3.2 Harmonic moments as commuting flows 



Equation (|3.5|) suggests to treat log/ a (;z) as a generating function of commuting flows 



with respect to spectral parameter a. The expansion 

hg f a (z) =H (z) + Y l (Ma)Hk(z)-MajHk(zj) (3.8) 



fc>i 



defines generators H k , H k of the commuting flows. Clearly, H (z) = log f(z). It implies 
evolution equations for f(z), 

dlogf(z) dH k (z) dlogf(z) dH k (z) 



dt k dt dt k dt 

and integrability conditions: 

dH k (z) _ dHj(z) dH k (z) _ dHj(z) dH k (z) dHj(z 



dtj dt k dtj dt k dtj dt k 



(3.9) 



The real part of ( |3.8|) vanishes on the boundary (as it is the Dirichlet Green function), 
therefore, the boundary values of H k and H k are complex conjugated: 



H k (z) = H k (z), ze-y. (3.10) 

The structure of integrable hierarchy becomes explicit if instead of functions of z one 
passes to functions of its image w under the map f Zo : w = f Zo (z) = f(z). 

Using the chain rule, one can write 

V(a)log/ 6 (z) = V(a)logf b (z(w)) + (V (a) log f(z))wd w log f b . 

w 

In the last term we observe that V(a) log /(z) = <9 io log f a (z) (using ( p.5|) at b = z ). 
Subtracting the same equality with a, b interchanged, we come to equation of the zero- 
curvature type: 

V(a) log f b - V(6) log f a - {log f a , log M = 0, 

where the Poisson brackets are defined as {/, g} = w*^—^- — w^-^- and t^-derivatives 

aw ot ow ot 

are taken now at fixed w. 

Let z(w) be the map inverse to w = f(z). Equation ( |3.5| ) at b = zq, being rewritten 
for the inverse map, has the form of a one-parametric family of evolution equations of 
the Lax type labeled by the spectral parameter a. They are 

V(a)z(w) = {logf a (z(w)) 7 z(w)}. (3.11) 

We refer to them as to deformation equations. The zero-curvature conditions ensure 
that these equations are consistent, i.e. flows with different values of spectral parameter 
commute. 
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The integrability conditions ( |3.9|) in the new variable acquire the form of the zero- 
curvature equations 

d tj Hi(w) - d ti H 3 {w) + {Hi(w), H 3 (w)} = 0, 

di.Hiiw) - d- u H 3 {w) + {Hj(w), Hi(w)} = 0, (3.12) 

d tj H t (w) + d^Hjiw) + {Hi(w), Hj(w)} = 0. 

From ( |3.6| ) it is easy to see that Hk are polynomials in w while Hk are polynomials in 
w^ 1 . Furthermore, for any basis such that ipk( z ) = 0( (z — Zo) k ) these polynomials are 
of degree k. In other words, the generators are meromorphic functions on the Riemann 
sphere with two marked points at w = and w = oo. This is a particular case of the 
universal Whitham hierarchy |[2| , known as dispersionless Toda lattice. In it is 



proved that the full set of zero- curvature conditions (P-12Q , together with the polynomial 
structure of the generators, already imply existence of the Lax function and Lax-Sato 
equations. 

For a particular choice of the basis, the Lax function can be expressed through the 
inverse conformal map. Consider the interior problem, set Zq = G D and fix the basis 
of holomorphic functions in D to be the natural one: ipk{ z ) = z k /k. From ( |3.7| ) and 
( |3.8|) it follows that Hk are holomorphic in D while are meromorphic with the fc-th 
order pole at so that Hk = z~ k + 0(1) as z — > 0. Combining these properties with the 
polynomial structure of the generators as functions of w, and taking into account (|3.10| ), 
one gets 

H k =(z-\w)) +Uz-\w)) 



(3.13) 

Hk=(z-\ w -^) < y-(z-\ w -^) o , 



where z(w) is the map inverse to f(z) and z(w) = z(w) is its Schwarz double. The sym- 
bols (/(iy))>o, (/(w))<o mean truncated Laurent series, where only terms with strictly 
positive (negative) powers of w are kept, while (/(w))o is the free term of the 
series. The free terms in ( |3.13|) are found with the help of an easily proved identity 

log = ^2 y( z ~ k ( w ))o zk valid for any series of the form w(z) = z + W\Z 2 + . . . . 

z k>i k 

Similar formulas hold for the exterior problem with the choice zq = oo, ipk{z) = z~ k /k. 



3.3 Deformation equations and dispersionless Toda hierarchy 

Expending the deformation equations ( |3.11|) in spectral parameter we obtain the Lax 
representation of the dToda hierarchy. To see this, we recall that the dToda hierarchy is 
an infinite set of evolution (Lax-Sato) equations for two Lax functions 



L(w) = rw + UkW 



-h 



k>0 



L(w) = rw 1 + ^2 UkW k 

k>0 
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The equations are 

dL(w) dL(w) ~ 

= {H k {w), L{w)} , —jfi— = i L M, H k {w)} 

and the same for L(w). Here the Poisson brackets with respect to w and to are defined 
as above, and 

H k (w) = (L k (w)) >o + 1 -(L k (w)) o 

H k {w) = (L\w)) <Q + l -(L\w)) o 

are generators of the flows. One may also introduce Hq(w) = log w. They obey the disper- 
sionless zero-curvature conditions ( |3.12| ) which express consistency of the Lax equations 
and generate an infinite set of partial differential equations for coefficients of the Lax 
functions. 

To summarize, we have the following identification of the Lax functions with conformal 
maps: 

• For the interior problem (z = with the natural basis): 

L(w) = ^— LH = -^r, (3.14) 

where z{w) is the function inverse to w = fo(z); 

• For the exterior problem (zq = oo with the natural basis): 

L(w) = z{w) , L(w) = ziw" 1 ) , (3.15) 
where z{w) is the function inverse to w = fooi z ) 

So the inverse map, z(w), and its Schwarz double, z^w^ 1 ), both obey the Lax equations. 

Let us comment on another choice of basis. Suppose A -1 is a local parameter at the 
normalization point zq, i.e., 

A = — h c + c^z-Zq) + 0( (z-z ) 2 ) , z ^ oo 

Z-Zq 

A = z + c + ^- + 0(z~ 2 ) , z = oo. 

with some domain- independent coefficients Cj. Assuming that A -1 (z) is a well defined 
local parameter in a domain that contains D, we choose the basis of holomorphic functions 
in D to be ipk{z) = \~ k (z)/k. This results in a linear change of times with the help of a 
triangular matrix. Repeating the above arguments, one identifies the Lax function with 
X(z): L(w) = X(z(w)). This Lax function provides a solution to an equivalent hierarchy 
in the sense of ||14|| . 
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3.4 String equations 



A customary way to fix a solution to a dispersionless hierarchy is to impose an additional 



constraint on the Lax functions (sometimes called string equation [13). We are going to 



show that not only the deformation equations but the string equation, too, is an easy 



consequence of the Hadamard formula and our basic relation (3.5). 

By r = r ZQ denote the conformal radius of the curve 7 with respect to the normaliza- 
tion point z . Let us calculate 5 e ^ logr in two different ways. The first one is to use the 



Hadamard formula (3J.) when both arguments tend to the normalization point: 



6e(ologr = ^\dsG(z ,0\ 2 = ^m(0\ 2 

(recall that |/(£)| = 1 for £ € 7). The second one is obtained from ( |3.3| ) in the limit 

zi^(£7, Z2, z 3 -> z : 

S e{0 logr = --d h log |/(0I, 
where we have used ( |2.3| ). Combining the results, we obtain the relation 

d to log\f(z)\ 2 = -\d z f(z)\\ ze 1 . 

Passing to the variable w = f(z), one rewrites this equation as 



^ I dz(w) dz(w)\ _ , . 

2Ke ( m ^^r) = 1 ' w = 1 ' 

where z(w) is the inverse map, as before. Being analytically continued from the unit 
circle, it reads 

{z(w), ziw- 1 )} = 1. 

This is the customary form of the semiclassical string equation. Similar arguments in 
case of the exterior problem lead to the same equation. 

For the natural basis, let us rewrite the string equation in terms of the Lax functions 
L,L. Using the identifications (|3.14j ) and ([3.15|) , we have: 



{L 1 (w), L 1 (u^)} = —1 for the interior problem 

(3.16) 

{L(w), L(w)} = 1 for the exterior problem 

So, although the string equation in terms of z(w) is the same, the interior and exterior 
problems correspond to two different solutions of the same (dToda) integrable hierarchy 
(cf.@). 



4 Tau-function and dispersionless Hirota equations 
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4.1 Tau- function 



Symmetry relation ( [3. 3D implies that there exists a real-valued function of harmonic 
moments F(t , t, t) = F(t ; t\, t 2 , . . . ; £1, i 2 , ■ ■ ■) such that 



G{zi, z 2 ) = g {z h z 2 ) + -V(z 1 )V(z 2 )F : 



where the function g does not depend on the domain, i.e., on moments (the coefficient \ 
is set for future convenience). Note that adding to F a quadratic form in £/. amounts 
to changing the function g only. 

From the definition of the Green function it follows that go{z\, z 2 ) is symmetric and 
harmonic in D in both arguments with the only singularity log \ z\ — z 2 \ as Z\ — > z 2 , or, in 
terms of a local parameter A -1 , log |A _1 (zi) — A - 1 (^2) | - In case of the interior problem, D 
can be any domain not containing the point of infinity. The function go must be the same 
for all such domains. Therefore, singularities of go in the extended complex plane, other 
than the logarithmic singularity at merging points, may occur only at infinity. In other 
words, g is of the form go(zi,z 2 ) = log |A -1 (zi) — A _1 (z 2 )| +h (z 1 , z 2 ) where h is regular 
and harmonic in both arguments everywhere in the complex plane but at infinity. The 
function ho is in fact a matter of definition of the tau-function. For some particular basis 
we are free to choose it to be zero. Hence we obtain the following important relation: 

G(zx, z 2 ) = log \X~ 1 (z 1 ) - X~ 1 (z 2 )\ + ^V(zi)V(z 2 )F , z h z 2 e D (4.1) 

Changing the normalization point or passing to another local parameter results in mod- 
ifying the function g . As is already pointed out, this is equivalent to adding to F a 
quadratic form in times. Taking into account the remark at the end of Sec. 3.4, one 
sees that this agrees with the relation between tau-functions of equivalent hierarchies 



discussed in JL7]] . It is easy to see that in this form the relation holds true for the exterior 
problem as well. 

For clarity, let us specify the above equation for the natural basis. Instead of the dif- 
ferential operator V(z) defined for arbitrary basis, it is convenient to use its specialization 
for the natural basis 

Viz) = d to + £ (^-d tk + ^dt k ) (4.2) 

so that V(z) = T>{z) for the exterior problem and V(z) = P(z _1 ) for the interior one. 
Eq. ( |4.1| ) is then rewritten in the form 

log \z\ — z 2 \ + -V(z^ 1 )V(z2 1 )F for the interior problem 
G(z 1 ,z 2 ) = { (4.3) 
log \zi 1 — z^ x \ + -T>(zi)V(z 2 )F for the exterior problem 

The function F plays a central role in what follows. It is the tau-function of curves 
discussed in [§, |3|, ||. It is a dispersionless limit of logarithm of tau-function of an 



integrable hierarchy or dispersionless tau-function |]12[ . For brevity we refer to it simply 
as tau-function. A relation between the tau-function and the Green function of the type 
(4.3) has been emphasized by L.Takhtajan In sect. |5] we give a few explicit integral 



representations for the tau-function. 
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4.2 Conformal maps and tau-function 

In this section, we use the natural basis. To unify formulas for the interior and exterior 
cases, we write ipk{z) = X~ k (z)/k, where X(z) = z~ l for the interior problem and X(z) = z 
for the exterior one. (Correspondingly, zq is either or oo.) We also employ the notation 

f(zj) = fj, X(zj) = Xj, 

so that V{z) = d to + D(z) + D(z). 

Representations of the conformal map through the tau-function can be obtained from 
( [4.3|) by separating holomorphic and antiholomorphic parts in z\ and Z2 of the Green 
function. Holomorphic parts in both variables yield 



h-h 

X\ — A2 2 



log = D(X 1 )D(X 2 )F - -d\F. (4.5) 



Terms holomorphic in z\ and antiholomorphic in Z2 yield 

log - jL^j = -D(X l )D(X 2 )F. (4.6) 

Various specifications of ( |4.3| ), ( [4.5[ ) and ( [4.6| ) lead to representations of the conformal 
maps through the tau function. Tending, for example, z 2 — > z in ( |4.3|) and extracting 
holomorphic parts in Z\ , we get the following expressions for the conformal maps through 
the tau-functions: 

f(z) = X(z) exp (- \d\F - d t0 D(X(z))F). (4.7) 
The leading coefficients as z — > Zq give formulas for the conformal radius: 

21ogr = -«ag ) F (4.8) 

(recall that the sign factor k is introduced in Sec. 2 to distinguish between interior and 
exterior problems). Limits z 2 — > z Q in (|4.5|) , ( |4.6|) give other representations for the 
conformal maps, 

f(z) = e-^° F (X(z) - dl t F - D(X(z))d tl F), 

(4.9) 

= e -^o F D{X{z))d h F 



and another representation of the conformal radius: 



r- 2K = d 2 tl - t F (4.10) 



These formulas agree with ( fD3| ) provided F satisfies the dispersionless Toda equation 

dl- t F = e d l F . (4.11) 



12 



In fact all coefficients of the Taylor expansion of the conformal map around the normal- 
ization point can be easily read from (fO|). 

Merging the points Z\ and Z2 in (|4.5|) , (|4.6|) , we get formulas for the derivative and 
the modules of the conformal maps: 

f'(z) = \\z)e D2{x{z))F -^ F 

(4-12) 

l-\f(z)r = e W (-D(\(z))D(\(z))F) . 



It is easy to see that the generators of commuting flows introduced in ( |3.8|) are ex- 
pressed through the tau-function as follows: 



H k {z) = \\z)-(±d tQ +D{\{z)))0, l l\ 



H k (z) = (±d t0 + D(\(z)))dz k F, 
H (z) = log A(*) - (U + D(X(z)))d t0 F. 



2 

As is seen from ( |4.3|) , coefficients of the Laurent series of the regular part of the Green 
function are second order derivatives of the tau-function. The formula 

2feflw+*4.)G(*,*)- 1 1 dl °^M dl °^M 



2ixi JdD d log z 



can be easily proved by substituting (|3.4j ) and evaluating residues at the poles 2:1,22- 
The left hand side of this formula could be thought of as a dispersionless analog of the 
Gelfand-Dikii resolvent. Expanding both sides in power series in Z\,Zi using ( }4.3|) and 
Q3.8|), we obtain 



d 2 F 1 r dH n dH m 



dt n dt m 2ni(n+m) 7 7 dlogz 



d 2 F 1 r dH n dH m 



n > 0, m > 1 



n, m > 1 



dt n dt m 2mn T~ t dlogz 
(in the last formula we used ( |3.10|) ). 

In addition to the relations above, we point out that the tau-function provides new 
representations for some classical objects of complex analysis. Consider, for example, the 

Schwarz derivative S(f, z) = -jrhj — § \ f'(z) J °^ ^ ne con f° rma l map /. A non-singular 
part of lim d z d^G(z, () is j^S(f, z), so equation ( [L3l) , with z = 00 and X(z) = z, yields || 

S(f,z) = Q(d z D(z)) 2 F 
where d z D(z) = z ~ l ®t k - The reproducing kernel fll| of the (exterior) domain D, 

k>l 

K(z, C) = 2d z d{G(z, C), is given by 

K{z,£) = d s D(z)d<D(t)F. (4.13) 
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4.3 Dispersionless Hirota equations 



The tau-function obeys an infinite set of non-linear differential equations. They can be 
obtained by excluding / from the identities of the type ( |4.5| , |4.6| ) by means of eq. (|4.7|) . 
Equations obtained in this way are dispersionless limits of the Hirota bilinear relations 
fl~8|| obeyed by the tau-functions of integrable hierarchies. 



For example, consider eq. (|4.5|) for pairs of points (zi,z 2 ), (z 2 ,z 3 ) and (z 3 , Zi). Expo- 
nentiating and summing the equations one gets the identity 



D{z 2 )D{z 3 )F 



D(z 3 )D( Z1 )F 



0. 



(4.14) 



obeyed by the tau-function F for both interior and exterior problems. This is the sym- 
metric form of Hirota's equation for the dKP hierarchy. It encodes algebraic relations 
between the second order derivatives of the tau-function. These relations are obtained 
on expanding the Hirota equation in powers of Zi,z 2 , z 3 and comparing coefficients^ The 
operators D(z) and D(z) are defined in (|4.4j) . Equation (|4.14|) , as well as other Hirota 
equations given below, are the same for interior and exterior problems. 

More general equations obtained in a similar way include derivatives with respect to 
to and ik- These are equations of the dToda hierarchy: 



Izi - z 2 e 



D(z 1 )D(z 2 )F 



-D{ Zl )D{z 2 )F 



Zie -9 t0 D( Zl )F _ Z2e ~O t0 D(z 2 )F^ 
J_ p dt (d t0 +D(z 1 )+D(z 2 ))F 



Z\Z 2 



(4.15) 
(4.16) 



Note that eq. ( |4.15| ) can be obtained from ( 4.14 ) in the formal limit z 3 — > if to un- 
derstand \im z ^ D(z) as — d to . These equations allow one to express the second order 
derivatives d? , F, d? r F with m,n > 1 as certain functions of the derivatives dl f F, 
dl o i k F. The dispersionless Toda equation ftlll|) is the limit of (|4.16|) as z\,z 2 — > oo. Each 
side of this equation is the squared conformal radius as one can see from ( |4.8| , |4.10|) . 

In the same manner one can derive other equivalent forms of the Hirota equations: 



ziiz?, - z 2 )e 



-Vizs^Diz^-Dizs^F 



Z2(.z 3 - zi)e 

z J Zl _ Z2 y(D( Z1 )-D(z A )){D(z 2 )~D(z z ))F 
1 _ e -(D(zi)-D(z 3 ))(D(z 2 )-B(z 3 ))F = 



V(z 3 )(D(z 2 )-D(z 3 ))F 



— {z l Z 3 ){z 2 ^-3 

where the operator T>(z) is defined in ( [4.2| ). Eqs. ( |4.17| ), ( [4.18| ) may be interpreted as 
analogs of eqs. ( |4.15| , [4.16| ) with the normalization point moved to z 3 , in which case 
d to should be substituted by V(z 3 ). In the limit z 3 — > oo they convert into eqs. (|4.15|) , 
( 4.16| ) respectively. In ( 4 . 1 8|) z 3 can be regarded as an independent formal variable, not 
necessarily complex conjugate to z 3 . We stress that the full set of differential equations 
for F obtained by expanding the Hirota equations ( |4.15| , [4.16|) or ( |4.17|J4.18| ) is the same. 

The special case of ( 4.1 8| ) z 2 = z\ is worth mentioning: 



z - l\ e V(z 3 )(V(z 3 )+D( Zl )-D(z 3 )+D(z 2 )-D(z 3 ))F 



(4.17) 



(4.18) 



-(D(z)-D(0)(5(z)-D(Q)F 



\z~ l -C 



-l\2 e V(QV{z)F _ 



2 In (T^, poj , the Hirota eq uatio n for the dKP hierarchy was obtained in an equivalent but less sym- 
metric form. It follows from (4.14) in the limit z% — » oo. 



14 



Further specialization £ — > z yields 

\z\%D(z)d,D(z)F = e v2{z)F , (4.19) 

This equation is especially remarkable. On the one hand, it looks like the dispersion- 
less Toda equation ( [4.1 1[ ) where d to , d tl , % are replaced by T>(z), z 2 d z D(z), z 2 dzD(z) 
respectively (moreover, ( |4.19|) becomes (|4.11|) as z — ■> oo). On the other hand, Q4.19Q is 
equivalent to the Liouville equation. Indeed, eq. ( [4. 19|) tells us that the field 



X (z)=V 2 (z)F-\og\z\ 4 

obeys the Liouville equation 

d t dgx = 2e x 

for z G D. (Here we imply that D contains oo, but similar equations can be written for 
the interior problem, too.) By virtue of ( |4.3| ), ( |4.12|) , the solution to this equation can 
be written as follows: 

X (z)=2hm [G{z,z')-\og\z-z'\\ 

z'^z 

or 

\f'{z)\ 2 



(l/WI 2 -!) 2 

Note that ds 2 = e x dzdz is the pull back of the Poincare metric with constant negative 
curvature = —2e~ x d z d z x — ~ 4 in D. We also note that the Liouville field x 

equals the value of the reproducing kernel |llj of the domain D at merging points: 



x(z) = \ogK(z, z), where K(z,() is the reproducing kernel ( |4.13| ). 



4.4 Residue formulas 

In this section we present formulas for third order derivatives of F. Formulas of this type 
are known in the theory of dispersionless integrable hierarchies and are referred to as 
residue formulas pi. They are used, in particular, to prove the associativity equations 



for tau-functions of Whitham hierarchies [21. 22. 23]. 



The basic relation to derive the residue formulas is (|3.2| ). By virtue of ( |4.1|) , its left 
hand side is a generating function for third order derivatives of the tau-function. Let us 
first rewrite this formula in holomorphic terms. To do that we note that d n G(z,^) = 
—2\d^G(z,^)\ for £ on the boundary (the vector gradG is normal to the boundary). 
Further, since G vanishes on the boundary, d^G(z,^,)d^ + d^G(z,^,)d^ = 0. Therefore, 
\d^Gd^\ = —nid^Gdl; holds on the boundary. Taking all this into account, we rewrite 
©as 

V(a)V(6)V(c)F = A / d t G(a,z)d z G(b z)d z G(cz) ( ^ )3 
ni j 7 dzdz 

Finally, using G(a,z) = — log |/ a (z)|, where f a (z) is a holomorphic functions of z (see 
(HP), we get 

V(aW(bW(c)F = - — I d\o&fa{z)d\ogf b {z)d\ogUz) 
y ' w y ' 2m L dzdz 
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Let us expand the both sides in the basis ipk i n eac h argument using (|2.2|) , (|3.8| ). Com- 
paring the coefficients, we come to the residue formulas for third order derivatives of the 
tau-function: 



d 3 F 



dtidt m dt n 
d 3 F 



1 f dHidH m dH n 
2%i J7 



dzdz 



1 

2ni j-y 



dHidH m dH n 



l,m,n > 0, 



L m > , n > 1. 



dtidt m dt n 2ni J-y dzdz 
The formulas ( f4 . 2 0| ) were first obtained by I.Krichever 15) by other means. 



(4.20) 



5 Integral representations of the tau-function 

In this section we discuss integral representations of the tau-functions and connection 
between the tau-functions for the interior and exterior Dirichlet problems. We employ 
the special notation D int for the interior domain bounded by the curve 7 and D ext for 
the exterior one. We use the natural basis, i.e., the normalization points are G D[ nt 
and 00 G D cxt , and the basis functions are z k /k and z~ k /k respectively. At last, the 
tau-functions are denoted as F mt , F ext . 



5.1 Electrostatic potentials 



Let us begin with the exterior problem. Our starting point is eq. ( |4.3| ) which we rewrite 
here introducing the auxiliary function (potential) 



ext 



Eq. ( fL3| ) then reads 



$ ext (^) =V(z)F 



^)$ cxt (0 = 2G(2,0-21og 



(5.1) 



-1 



r 1 



for both z and £ being in D ext . (We recall that T>(z) = d to + D(z) + D{z).) This function 
admits an electrostatic interpretation explained below. When one of the points reaches 
the boundary, the Green function vanishes and one has 



ext / 



^logl^-r 1 !, £e 7 - 



The last formula makes sense of the variation of the function $ cxt under the bump 
deformation of the domain (see sect. 0). According to ( |2.3| ), the variation is 



-logi^-r 1 ! 

7T 



and, therefore, one may write the integral representation 



/ Alog^-C 1 



(5.2) 
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up to a harmonic function that do not depend on shape of the domain. Adding to (|5.2| ) 
harmonic functions independent of times amounts to redefinition of the tau-function 
by terms linear in times, which contribute neither to the Green functions nor to the 
formulas for conformal maps ( [4.7] ), ( f4.9p , and leave Hirota equations unchanged. We set 
all such terms to vanish. In a similar way, for the interior problem we get S e ^^ int (z) = 

— log \z — £| and 

7T 

^\z) = - [ d\\og\z-C\-c, 
vr Jfl cxt 

where the constant c is necessary for regularization of the divergent integral over D ext . Let 

2 f 

us cut off the integral at a big circle of radius R and set c = c(R) = — d 2 (\og \ C\ — 

7T J\£\<R 

R 2 log R 2 - R 2 . Then we get 



$ int (^) = lim (- [ d 2 C log \z- C|- c(R)) =- I efClogll-zC 1 !-^, (5.3) 

R^co VtT JD^t / 71 JD cx t 

where we introduced 



v = — [ log \z\ 2 d 2 z. 

7T J Dux 

Formulas (|5.2|) , (|5.3| ) present the potentials <3> ext , $ int in the form of an integral over 
the complementary domain, since by definition (|5.1|) the functions $ ext , $ mt are harmonic 
in -D ext , Ant respectively. One may use these formulas to extend the functions $ to the 
whole complex plane. For example, (5.2) defines a function $ cxt which is harmonic in the 
exterior domain and satisfies the Poisson equation — d z dz§ ext (z) = 1 — Tit^i^z) in D int . 
This function is the potential generated by a uniformly distributed charge in D int and a 
compensating point-like charge at the origin. The expansions of $ ext at small and large 
z read: 

$ cxt (^) = -\ z \ 2 + 2t log \z\ + ^(t^z* + iTz k ) , z -> 0, 

fc>0 

(5.4) 

$ cxt (z) = Vo + ^(4 nt z~ fc + ^k^~ k ) , z -> oo, 

k>0 

where the additional superscripts are set to distinguish between exterior and interior 
moments. Under our assumptions, these series converge in D int and D ext respectively. 
The functions $ cxt and d z $> ext are continuous at the boundary. 



5.2 Integral formulas for the tau-function 



Using the same strategy, we set z in ( |5.1| ) to the boundary and interpret this formula as 
a result of the bump deformation of the domain. It is easy to check that the variation of 



1 1 



IS 



F cxt = !-[ d 2 z^\z) = --J d 2 z [ d 2 ( log 

27T J D int 7l Z J D int J Dint 

-$ cxt (0 --J d 2 z\og\z- x -C x \ = -$ cxt (0- 

7l Z J D^t 7T 



e 

2^ 



(5.5) 
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Eq. C |5.5|) presents the tau-function for the exterior Dirichlet problem as a double integral 
over the domain complement to the D ex± [||]. 

Similar arguments give the tau-function of the interior problem: 

F int = lim (- -L / d 2 z [ d 2 (\og \z-C\ + C(R) - c(R)t ) 
with c(R) as in (|5.3|), and 

C{R) = -J d 2 z [ d 2 C\og\z-C\ = -R\2\ogR 2 -l). 

7T Z J\z\<R J\t\<R 4 

It is implied that the integral over D ext is cut off at R. Note the relation 

F iat + F ext = J_ f d 2 z f ^l 0g |l_^-l|2 ( 5 . 6 ) 



7T 2 JD int JD. 



'ext. 



In fact, the tau- functions admit other useful integral representations. Let us mention 
some of them. In terms of the potential $ ext , extended to the whole complex plane 
as explained above, the tau-function of the exterior problem may be represented as a 
(regularized) integral over the whole complex plane: 



^ cxt = lim ( 



e^O \ 2-K JC\D 



dJSf 



ext. 



d 2 z + tl\oge\ (5.7) 



where D £ is a small disk of radius e centered at the origin. In the electrostatic interpreta- 
tion, the tau-function is basically the energy of the system of charges mentioned above. 
Stokes formula gives an integral representation through a contour integral: 

™ vt 1 9 1 f ^ / s zdz — zdz 
Similar formulas can be written for the interior problem. 



5.3 Legendre transform F mt <->■ F cxt 

The tau-functions for the interior and exterior problems are connected by a Legendre 
transform. The former is the function of the interior moments £™* and the area to, while 
the latter is the function of the exterior moments t| xt and the area. From ( |5.1| ) we read 
that first order derivatives of the tau-function for the Dirichlet problem in D int or Z} ext 
with respect to the harmonic moments are (up to the factor k) harmonic moments of the 
complimentary domain: 

f) pint a rpext 

ktT\ -5=5- = ^ (5-9) 



and 



<94 nt ' <9t e fc xt 



= V . 



dt dt 
Using (|5.6| ), we obtain the relation 



F int (t , t int , t int ) = J2 ( k 4 nt 4 xt + k 4 nt - ^ cxt (^o, t ext , t 



ext\ 



k=l 
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By virtue of ( |5.9|) , it means that the functions F mt and F ext are related by Legendre 
transforms with respect to all variables but to'- 



Mat = £ t ext ^_ £ -ext ^_ _ ^ext 

fc=i V ol k ol k ) 



°° ( f)F int flF int \ 

pent _ \ " / Ant u± _|_ rint ul 1 _ pint 



(5.10) 



5.4 Homogeneity properties of the tau-function 

Expansion of the potential $ ext (;z) around the origin ( |5.4j ) allows one to prove another 
important property of the tau-functions. Integrating both sides of it over D int and using 
( |5.9|) , we obtain a quasihomogeneity condition for _F ext . A similar condition for F mt is 
most easily derived from ( |5.10| ). They are: 

4F int = t 2 + 2t d to F mt + E(2 + k){tkdt k F mt + 4^F int ), 

(5.11) 

fc>0 

These formulas reflect the scaling of moments as z — > A2 with real A: tf* -> X 2+k tf\ 
tf^ — > A 2 ~ fc t 1 fc nt . The logarithmic moment t>o, under the same rescaling, exhibits a more 
complicated behaviour: v — > A 2 f + to-^ 2 log A 2 . To get rid of the "anomaly term" £q one 
may modify the tau-function by subtracting ^logf 2 ,. 



6 Dirichlet problem on the plane with a gap 

Consider the case when the domain D mt shrinks to a segment of the real axis. Then the 
interior Dirichlet problem does not seem to make sense anymore but the exterior one 
is still well-posed: find a bounded harmonic function in the complex plane such that it 
equals a given function on the segment. The problem admits an explicit solution (see e.g. 
P^fl). Possible variations of the data are variation of the function on the segment and 
the endpoints of the segment. Solution of this problem as well as its integrable structure 
may be obtained from the formulas for a smooth domain as a result of a singular limit 
when a smooth domain shrinks to the segment. 

The tau-function, obtained in this way, is a partition function of the Hermitean one- 
matrix model for the one-cut solution in the planar large iV limit fH 



6.1 Shrinking the domain: the limiting procedure 

Let us consider a family of curves 7(e) obtained from a given curve 7 by rescaling of the 
y-axis as y — > ey. (If the curve 7 is given by an equation P(x, y) = 0, then 7(e) is given 
by P(x,y/e) = 0.) We are interested in the limit e — > 0, 7(0) being a segment of the 
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V 



p(x) 

a (3 x 

Figure 2: A thin domain stretched along the real axis with thickness Ay(x) shrinks into 
a segment with density p(x) . 

real axis. We denote the endpoints by a, (3. Let Ay(x) be the thickness of the domain 
bounded by the curve 7(e) at the point x (see Fig. |]). 

We introduce the function 

p(x) = lim — ^ 
e->0 s 

It is easy to see that in case of general position this function can be represented as 

p(x) = J(x - ct)(0 - x) M(x), (6.1) 
where M(x) is a smooth function regular at the edges. So, instead of the space of contours 



7 we have the space of real positive functions p(x) of the form (|6.1| ) with a finite support 
[at, 0\ (the endpoints of which are not fixed but may vary). 

Let us use the first equation in ( |5.4|) to define times as coefficients of the expansion 
of the potential $(x) = $ ext (x) generated by the thin domain with a uniform charge 
density (and with a point-like charge at the origin). In the leading order in e the rescaled 
potential is 



x 



l im ^M = _i f dx'p{x')\og(---) 2 = T \ogx 2 + Y. T ^ k - (6-2) 

e^O £ 7T J a ^X X 1 J ^[ 



Comparing this with (|5.4| ), we get 

T = lime _1 t , 



T k = Jime- 1 ^ + if) , k > 1, k ? 2, (6.3) 

T 2 = lime" 1 ^ + ^ xt - 1). 
1 [P 

Note that Tq = — p{x)dx but similar integral representations for other times, = 

7T J a 

/ p(x)x dx, which formally follow from ( |6.2| ), are ill-defined. On the other hand, 

Ja 
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2 fP 
ilk 



2\ 



harmonic moments of the interior behave, in the scaling limit, as t k nt = ek l pk + 0(e 
where pk are well-defined moments of the function p on the segment: 

1 rP 

= — p(x)x dx . (6.4) 

7T Ja 



Using integral formula ( |5.5| ), it is now straightforward to find the scaling limit of the 
tau-function F ext for the exterior problem. Taking into account ( |6.3|) , we introduce the 
function F cut as follows: 

F cxt (et ; eti, \+et 2) et^, . . . ; e?i, ^+ei 2 , sts, . . . ) = 

(6.5) 

= e*F™\T Q] T u T 2 , T 3 ,...) +0(e 3 ), 

where To = to, T& = £& + ife. Since the second order derivatives of F ext are invariant under 
the rescaling ( |6.5| ), the function T cut obeys the Hirota equations ( |4.14j) -( 4.16|) (as F ext 
does) where one has to set d tk = di k = dr k - The latter means that _F cut is a solution of 
the reduced dToda hierarchy (see e.g. [P~2]). This sort of reduction is usually referred to 
as dispersionless Toda chain. We conjecture that other types of reduction correspond, in 
the same way, to shrinking of _D int to slit domains of a more complicated form. 

An integral representation for the F cut (obtained as a limit of ( |5.5| )) reads: 

pent 

7T^ J a J a 



1 fP [P 

p(xi)p(x 2 ) log \x x — x 2 \dxxdx 2 

1 [P 1 

Represent this as F cut = — / p(x)<j)(x)dx = - V^/i^Tfe, where pk are defined in (|6~1 ) 

2tt J a 2 fe > 

and 

1 rP 

/A) = — / p(x) logx dx . 

7T Ja 

It is clear from the limit of (5^) that pu = dx k F cut for k > 0. Therefore, we obtain the 
relation 

2F cut = J2T k d Tk F mt } (6.6) 

fc>0 

which means that F cut is a homogeneous function of degree 2. This formula also means 
that the tau-function F cut is "self-dual" under the Legendre transform with respect to 

T , Ti, T 2 , However, the analog of the Legendre transform ( |5.10|) we discussed in 

sect. |y| does not include T , so the analog of the function F mt is the function 

npcut 

£ cut = ]TT fc/Ufc -F cut = F cut -T — — = 
k>i dI o 

rP rP 



1 fP ff> 

/ / p(xi)p(x 2 ) log \xi - x 2 \dx\dx 2 , 

7T Z Ja Ja 



which is the electrostatic energy of the segment with the charge density p on it, regarded 
as a function of the variables T , pi, p 2 , . . .. Properties of this function and its possible 
relation to the Hamburger ID moment problem are to be further investigated. 
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6.2 Conformal maps 

The conformal map / : C \ [at, 3] — > E from the plane with the gap onto the exterior E 
of the unit disk is given by the explicit formula 

2z - a - 3 + 2J(z-a)(z- 3) 

m l3 V a 

All formulas which connect conformal maps and tau-function remain true in this case, 
too. In particular, expanding f(z) as z — >• oo, 



read from (|4.9f) formulas for the endpoints of the cut: 

3 - a = 4 exp (\d 2 To F c ^) , 3 + a = 2d 2 ToTl F™\ 

The Green function is expressed through the f(z) by the same formula (|1.2j ). 
Set w = f(z), then the inverse map is 

3-a, _ x 8 + a 
ziw) = — - — ( w + w JH — . 

Clearly, z(u>) = so the constraint on the Lax functions is now L(w) = L(w) 

which signifies the reduction to the dispersionless Toda chain. 



6.3 Relation to matrix models 

In analogy with variations of the data of the Dirichlet problem discussed in Sec. 3, we 
can consider the following problem: given a set of Tk, k > 0, to find endpoints a, 3 of 
the segment and the density p(x) as functions of these parameters. This problem has 
ben appeared in studies of the planar iV — > oo limit of the Hermitean matrix model. In 



this case the function p{x) is a density of eigenvalues ||25|| . For completeness, we recall 
the basic points. 



Taking derivative of (|6.2|), we get the equation which connects the set of T k with a, 3 
and p{x): 

2 tP p(x>)dx> 

-v.p. / — = V(x), 

IT J a X ' — X 



where 



Consider the function 



V(x) = <f>(x) - T logx 2 = J2 T ^ k ■ (6.7) 



k>l 



W(z) = - f 

71 Jo 



I 3 p(x)dx 



z — X 



It is analytic in the complex plane with the cut [a,/3]. As z — > oo, it behaves as W(z) = 
T$z~ l +0(z~ 2 ). On the cut, its boundary values are W(x±i0) = —\V\x) =fip(x). The 
function W(z) is uniquely defined by these analytic properties. So, to find p one should 
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find the holomorphic function from its mean value on the cut. The result is given by the 
explicit formula 

4«r c-* v/(c- Q )k-/3) 

where the contour encircles the cut but not the point z. The endpoints of the cut are 
fixed by comparing the leading terms of (|6.8| ) with the required asymptote of the W(z) 
as z —>■ oo. This leads to the hodograph-like formulas 

1 i V'(z)dz _ Q 1 ^ zV'(z)dz = _ 9Tq 



2iriJ J( z - a )(z-P) 2m J ^{z-a){z-(3) 

which implicitly determine a, (3 as functions of T^. 

It follows from the above that F cnt coincides with the one-cut free energy of the 
Hermitean one- matrix model with potential ( |6.7| ) in the planar large N limit. (As a 
matter of fact, it was shown in |26|| that the partition function of the matrix model at 



finite N is the tau-function of the Toda chain hierarchy with dispersion.) Combining 
( |6.6| ) with the e — > limit of the second relation in (|5.11|) , we obtain 

J2kT k d Tk F cut + T* = 0. 

k>l 

This identity is known as the Virasoro L -constraint on the dispersionless tau-function 
F cut pB] , IT2J. Other Virasoro constraints can be obtained in the e — > limit from the 
W^i+oo-constraints on the tau-function _p ext . We do not discuss them here. 



6.4 An example: Gaussian matrix model 

If only first three variables are nonzero (i.e., T ,T , 1 ,T 2 ^ 0) function F cut can be found 
explicitly. Consider a family of ellipses with axes / and es centered at some point xq on 
real axis: 

{x - xq) 2 y 2 I 



I 2 e 2 s 2 4' 

The harmonic moments, as e — > 0, are (see Appendix in ref. [Q]): to = h^sl, t^ xt = 
2exosl~ 1 + 0(e 2 ), tf^ — 5 — esl" 1 + 0(e 2 ) and all other moments of the complement to 
ellipse vanish. From (|6.3|) we read values of the rescaled variables: 



T* = \{p-a)s, T 1 = 2^±^s } T 2 



4 (3 — a (3 — a 

and all the rest are zero. Here we have expressed times through the endpoints a = xo — ^l, 
(3 = x + |Z of the segment and the extra parameter s. The density function is 

p(x) = -T 2 ^J(f3 - x)(x - a). 

(Note that T 2 < 0.) 

Using explicit form of the tau-funtion for ellipse 0. pi . 



next 



1,2 1 ^0 3 2 . , ^1*1 + ^2 + ^1*2 

= -tn lOff = t n + t() = " 

ellipse 2 ° 6 1 - At 2 t 2 4 ° U 1 - 4t 2 *2 
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and the scaling procedure (|6.5|) , we find 



F cut (T , 7^2,0,0, 



r 2 iog 



-2T 2 



3 T2 r 2? 



4 



4T 2 



This is expression for the free energy of the Gaussian matrix model in the planar large 
N limit [p5fl. The same result can be obtained from the integral formula for F cut . 
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